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ABSTRACT 

Frequency dependence of pulsar linear polarization is investigated by simulations of 
emission and propagation processes. Linearly polarized waves are generated through 
curvature radiation by relativistic particles streaming along curved magnetic field lines, 
which have ordinary mode (O-mode) and extra-ordinary mode (X-mode) components. 
As emitted waves propagate outwards, two mode components are separated due to re¬ 
fraction of the O mode, and their polarization states are also modified. According 
to the radius to frequency mapping, low frequency emission is generated from higher 
magnetosphere, where significant rotation effect leads the X and O modes to be sepa¬ 
rated. Hence, the low frequency radiation has a large fraction of linear polarization. As 
the frequency increases, emission is generated from lower heights, where the rotation 
effect becomes weaker and the distribution regions of two modes are more overlapped. 
Hence, more significant depolarization appears for emission at higher frequencies. In 
addition, refraction effect of the O mode becomes serious in very deep magnetosphere, 
which bends the O mode emission towards outer parts of a pulsar beam and also causes 
the separation of mode distribution regions and hence the fractional linear polariza¬ 
tion increasing with frequency. If emission of different frequencies is generated from a 
region of the same height, serious O mode refraction can result in the decrease of both 
profile width and fractional linear polarization. The observed frequency dependence 
of linear polarization for some pulsars can be naturally explained within the scope of 
our scenario. 
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1 INTRODUCTION 


Pulsars are generally believed to be rapidly rotating neu¬ 
tron stars. Their polarized emission has been observed 
over a wide frequency range from about 100 MHz to 


32GHz over more than forty years (e.g. Manchester 

1971; 

Wu et al. 1993|: Xilouris et al. 1996|: Gould & Lvne 

1998|; 

Weltevrede & Johnston 2008; Han et al. 2009; Tiburzi 

et al. 

2013h. Very rich polarization features have been uncovered 


for both linear and circular polarizations, e.g., ‘S’-shaped 
position angle curves, orthogonal polarization modes, sin¬ 
gle sign and sign reversals of circular polarization, and their 
evolutions with frequency, etc. Linear polarization of aver¬ 
age pulse profiles is predominant for most pulsars at low 
frequencies. The polarization percentage gene rally dec reases 
with the increase of observing frequency (I Morris et al • 1970; 
iManchested Il97ll) . iManchester et, al] ll 1 97.T I noticed that 
some pulsars exhibit a critical frequency, below which the 
fractional linear polarization is constant, but above which 
it decreases with increasing frequency. To investigate the 
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chang e of pulsar polarization s across a much wider frequency 
range, IXilouris et al .11 (fl996 ) extended polarization observa¬ 
tions to mm-wavelengths, and notic ed sig nific ant d epolar- 
ization at high frequencies. Recently. I Johnston et al.l ((2008) 
pointed out that simple profiles are more likely to be char¬ 
acterized by high linear polarization than complex profiles, 
and polarization properties for different profile components 
of a given pulsar behave differently with frequency, e.g., PSR 
J0922+0638. 

In addition to analysis of polarization features for in¬ 
dividual pulsars, statistical investigations were also under¬ 
taken for pulsar linear polarization. It has been noticed at 
early days that s hort period pulsars tend to have high lin - 
ear polarization ([Huguenin et al.l Il97ll : I Morris et al.l Il98lh . 
Long period pulsars exhibit faster de polarization o ver fre ¬ 
quency than the sh ort period ones ([Morris et al.l Il98lh . 
IXilouris et al.l (|l995h presented the anti-correlation between 
the depolarization index and the acceleration potential 
near neutron star surface. Several authors also demon¬ 
strated that the degree of linear polarization is highe r for 
a pulsar with a large spin-down luminosity E ([Qiao et al.l 
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1995: von Hoensbroec h~et alJ Il998l : ICrawford et al.l l200ll : 
Weltevrede & Johnston 2008|). 


To understand various polarization features, numer¬ 
ous theoretical researches on pulsar polariz ation have been 
developed on th e emission processes (e.g. Xu et alJ I 2 OOO; 


Dvks et all l20ld : IWang et~ahl 120121 : iKumar &; Gangadharal 
2012 ) or the propagatio n effects (e.g. IWang et al.l l 20 ld : 
Beskin fc Philinnovl l2012h , both of which can lead to de¬ 


polarization. The depolarization of pulsar linear polar¬ 
ization is generally attributed to simultaneous interac¬ 
tion of two modes of orthogonally polarized radia tion 
([McKinnon &; Stinebring|[l998l : LKar astergiou et al.ll2005h or 
owing to randomization of position angles in weaker mag¬ 
netic fields in outer magnetosphere (Manches ter et al.lll975l : 
iMorris et al.l Il98lh . For example, McKhmon| (|l997l ) at¬ 
tributed the depolarization to birefringe nce of t he X mode 
and O mode above pulsar polar caps. Ivon Hoensbroech et ahl 
(119981 ) qualitatively explained depolarization at high fre¬ 
quencies by considering the evolution of one natural plasma 
mode within pulsar magnetosphere. However, the formation 
and distribution of the two modes within a pulsar beam were 
not mentioned; the refraction effect was simply discussed by 
analogy with the calcite crystal; the polarization limiting 
effect was simply analysed; and the evolution of linear po¬ 
larization across field line planes was not investigated. To 
get a better understanding of the frequency dependence of 
pulsar linear polarization, emission processes together with 
propagation effects need to be investigated systematically. 

Recently, we jointly studied the polarized curvature 
radiation together with propagation pr ocesses w ithin pul¬ 
sar magnetosphere numerically (IWang et al.ll 20 lJ j. We suc¬ 
ceeded in demonstrating the distributions of the X-mode and 
O-mode within a pulsar magnetosphere and explained the 
depolarization across an entire pulsar beam. We found that 
the depolarization is serious if the corotation of relativis¬ 
tic particles is not considered. When the rotation is taken 
into account, significant linear polarization will be produced. 
However, the frequency dependence behaviour of linear po¬ 
larization was not incorporated in the analysis. 

Stimulated by the observational and theoretical inves¬ 
tigations, we hereby develop the joint research of emission 
and propagation processes to understand frequency depen¬ 
dence of pulsar linear polarization by considering the polar¬ 
ized curvature radiation process together with propagation 
effects in a pulsar magnetosphere. In Section 2, we present 
theoretical basics for our calculations. Linearly polarized ra¬ 
diation within an entire pulsar beam, polarized pulse profile 
and their evolution with frequency are described in Section 

3. Comparisons with observations are elaborated in Section 

4. Discussions and conclusions are presented in Section 5. 


2 THEORETICAL BASICS 
2.1 Geometry and emission 

Pulsar radio emission is generally believed to be generated 
by relativistic particles streaming out along the curved mag¬ 
netic field lines within the open magnetosphere. Due to the 
bending of held lines, relativistic particles will experience 
perpendicular acceleration and produce curvature radiation. 
Coherent curvature radiation from bunches of relativistic 



Figure 1. Emission regions on a bundle of magnetic held lines 
with 0 ~ 2/3 6nm for a given sight line. The magnetic moment 
m is inclined by an angle of a with respect to the rotation axis 
f2. The thick line outlines all tangential emission points of held 
lines. 6nm represents the angle between the wave vector n and 
the magnetic momentum m. The grey area represents the region 
for outflowing relativistic particles along a bundle of magnetic 
held lines, s = 0.0 denotes the magnetic axis, s = 1.0 denotes the 
last open held line (LOF). 


particles serves as one of the most probable mechanisms for 
pulsar radio emission. The coherent particle bunch can be 
treated as a point-like huge charge for the sake of simplicity. 
Pulsar magnetic helds are assumed to be of the static dipole 
form in the radio emission region, 

B = R*( — ) 3 [3f (r • to) — to], (1) 

here B* is magnetic held strength at the magnetic equator 
of a neutron star, R* is the neutron star radius, r is the 
unit vector along r, and rh represents the unit vector of the 
magnetic dipole moment. Within the static dipole magne¬ 
tosphere, emission generated from the regions of the same 
polar angle 0 within a given magnetic held line plane will 
point towards the same direction, as shown in Fig. |T] 

Along a held line, emission radiated at different heights 
has different frequencies. The characteristic frequency for 
curvature radiation reads, 


Z'CR 


37 3 c 
47 Tp 


( 2 ) 


Here, 7 is the Lorentz factor of a relativistic particle, p repre¬ 
sents the curvature radius for instaneous particle trajectory. 
The larger curvature radius p for particle trajectory in the 
high magnetosphere naturally leads to lower frequency emis¬ 
sion. This is the physical basis for the radius to frequency 
mapping. 

For an inclined dipole magnetosphere, the pulsar emis¬ 
sion beam is generally compressed in the mer idional plan e 
dehned by the rotation axis and magnetic axis (|Biggj|l990l ) . 
For simplicity, our calculations assume that pulsar beams 
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and velocities of relativistic particles will be bent towards 
the rotation direction compared with those for the static 
dipole magnetosphere. Hence, the emission regions become 
asymmetry, with heights for the leading part of the beam 
getting larger than those for the trailing part, as shown in 
the right panels of Fig. [2j 

For emission of a single particle bunch at a given height 
in the pulsar magnetosphere, the radiation field E(t) and 
the corresponding Fourier components E(uj) can be calcu¬ 
lated b y using the circ ular path approximation, as described 
bv I Wang et al.l ( 2012 b The emission from a relativistic par¬ 
ticle bunch is beamed in a I /7 cone around the velocity 
direction. Therefore, the detectable emission from a given 
height and rotation phase is contributed by relativistic par¬ 
ticles not only at the tangential emission point of a field line 
but also on nearby field lines within the I /7 cone. By inte¬ 
grating the emission within an entire emission cone from all 
field lines, we will obtain the total polarized radiations in 
the given direction. In general, curvature radiation is highly 
linearly polarized, and the rotation h as great influences on 
polarized curvature radiation process (|Wang et al.ll20l3 b 



-5 0 5 -5 0 5 

H°) ri°) 

is = 3.0GHz without rotation is = 3.0GHz with rotation 


Figure 2. Emission heights (in the unit of R *) for different fre¬ 
quencies in a pulsar beam, plotted in the rotation phase ip and 
the sight line angle £. Corotation of relativistic particles is not 
considered for the left panels, but considered for the right panels. 
The parameters used for model calculations are an inclination an¬ 
gle of a = 30°, a Lorentz factor of 7 = 500, and a rotation period 
of a pulsar P = Is. 


are circular, i.e., the last open field lines which define the 
boundary of a pulsar beam have the sa me field line c onstant 
of r e ,i 0 f, as defined by equation (15) of iGangadharal (l2004l h 
In general, the beam size in the meridional direction does 
not affect our conclusions in this paper. 

The emission regions at three different frequencies of 
iscr are calculated as shown in Fig. [2] The height ranges are 
quasi symmetric about the magnetic axis for emission at a 
given frequency generated from the static dipole magneto¬ 
sphere. Emission near the beam centre comes from a lower 
altitude, as shown in the left panels of Fig. [2j because the 
curvature radii of the inner field lines at a lower altitude 
are comparable to those of the outer field lines at a higher 
altitude. Further, the entire emission regions get deeper in 
the magnetosphere for a higher frequency, which is caused 
by the smaller curvature radii of field lines in the deeper 
magnetosphere. It should be noted that relativistic particles 
within a pulsar magnetosphere not only stream along the 
curved magnetic field lines, but also corotate with magneto¬ 
sphere. When the corotation is considered, the trajectories 


2.2 Propagation 

Once polarized emission is generated through the curvature 
radiation process, they will be coupled to local plasma modes 
to propagate out in a pulsar magnetosphere, which is filled 
with relativistic plasmas that can be simply assumed to be 
cold (with a single Lorentz factor 7 ) and with a density 
of N p — t]Ngj- Here, 77 is the multipli city factor and TVg j 
repre sents the Goldreich-Julian density (|Goldreich &; Julianl 
[1969( 1. Two eigen transverse modes, the X-mode and O- 
mode, of the plas ma have been carefully i nvesti gated by 
I Wang et al.l (201 0) and lBeskin fe Philippov (12012 b The X- 
mode waves have a refraction index of nx = 1 and propagate 
rectilinearly; while the O-mode w aves have a refraction in - 
dex of no < 1 and suffer refraction ([Barnard fe Aronslll986h . 
The polarization states of the X-mode and O-mode change 
during their propagation. In the inner magnetosphere near 
the emission region, the polarization states of both modes 
evolve adiabatically, i.e., with the polarization of the X- 
mode wave orthogonal to the local k — B plane and the 
polarizat ion of the O -mode wave within the k — B plane 
(jCheng &; Rudermanll 19791 ). As the waves propagate to the 
higher magnetosphere, the adiabatic condition is not satis¬ 
fied and wave-mode coupling happens, i.e., one mode leaks 
to the other and vice versa. The altitude where the mode 
coupling happens is named as being the polarization limiting 
radius ([Cheng &; Rudermanlll979l : IWang et al Above 

the polarization limiting radius, the polarization states of 
both modes are frozen and will not be further affected by 
the plasma, except for possible cyclotron absorption in the 
outer magnetosphere. 

In addition to evolution of a single wave mode, the dis¬ 
tributions of polarized waves within each I /7 e missio n cone 
also c hange during propagation. As shown by I Wang et al.l 
(|2014h . the O mode refraction will lead detectable X and O 
modes at a given position to be incoherent, as they origi¬ 
nate from different parts of a pulsar beam. Adiabatic walk¬ 
ing will result in the ordering of polarized waves within the 
I /7 cone, i.e., with the O-mode wave vectors following the 
polarization plane of initial emission, but wave vectors for 
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Figure 3. Pulsar polarization emission beams at frequencies of v =0.8GHz, 1.0GHz, 1.4GHz, 2.0GHz, 2.4GHz and 3.0GHz from top to 
bottom for the X-mode intensities (/x), the O-mode intensities (/o), the total intensity (/x + /o), their difference (/x — Io, i.e., the 
net linear polarization) and the fractional linear polarization (L% — |(/x — Io)/(Ix + /o)| X 100%). The fractional linear polarization 
is plotted for the regions where the total intensity exceeds 1% of the peak intensity. The density distribution for relativistic particles 
is assumed to be located in a cone at s p = 0.5 and with a width of cr s = 0.12. The other parameters used for model calculations are 
a = 30°, 7 = 500, 77 = 100, B* = 10 12 Gs and P = Is. 


the X-mode rotate by 90°. Therefore, incoherent superposi¬ 
tion of the orthogonal X-mode and O-mode emission causes 
depolarization for total radiation power. If the corotation of 
relativistic particles is not considered, significant depolar¬ 
ization occurs. When the corotation is considered, the net 


polarization will reserve polarization features of the stronger 
mode. 
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(a) C = 30° (b) C = 31° (c) C = 32° (d) C = 33° 

Figure 4. Polarized pulse profiles resulting from cutting of the beams in Fig. [3] by the sight lines at f = 30°, 31°, 32° and 33°. The 
solid lines are for the total intensities, /x +^CO the dashed lines for the linear polarization intensities, |/x — /oh and the dotted lines for 
position angles of linear polarization. The polarization modes (X or O) are marked near position angle curves. 


3 EMISSION AT VARIOUS FREQUENCIES 

After knowing emission regions and propagation effects, we 
can now calculate the polarization states of emission at var¬ 
ious frequencies. The wave mode distributions within an en¬ 
tire pulsar beam (the X mode intensity, ix, the O mode 
intensity, To, total intensity, /x + /o, and the intensity for 
the linear polarization, /x — io), polarized pulse profiles, 
and their evolution with frequency can be calculated for the 
emission generated from different heights for the entire open 
magnetosphere. 


3.1 Emission generated by particles with a single 
value of the Lorentz factor 7 at uck 

In our calculations, the coherent particle bunch is simply 
treated as a point-like huge charge. The emission inten¬ 
sity is closely related to local particle densities. Here, the 
conal-shaped density distribution for particles is adopted as 
N(r,0,(f>) = N 0 r~ 3 f{9)g((f>), with 

fO) = fo exp[— 2 Jp) ^ 

5 ( 0 ) = 1 - ( 3 ) 

Here, s = 0 c /0 c ,max , 0 C is the polar angle of a field line footed 
on the neutron star surface, 6 c ,max = sin r e ,iof) rep¬ 

resents the polar angle maximum, i.e., the polar angle for the 
last open field line footed on the neutron star surface, s p is 
the position for density pealQ 

Using such a particle density distribution, we can now 
get the polarized emis sion at various frequenc ies in the 
model. The steps follow IWang et ahl (|2012l . l2014h . First, we 
find the field lines for emission within the I /7 cone point¬ 
ing towards a given direction. Then, the electric fields of 
the emission within the cone will be decomposed into the 
X-mode and O-mode components following the direction of 
local magnetic field. Evolution of the X-mode and O-mode 

1 In this paper the meaning of ‘s’ is the same as $ in I Wang et al.l 
(l2014l h 


components within the I /7 cone will be traced until the po¬ 
larization limiting radii by considering the refraction and 
the adiabatic walking effects. The integration of emission 
within each cone will result in the X-mode and O-mode in¬ 
tensities. Finally, emission pointing outwards from the entire 
open field line region will be mapped. The equations govern¬ 
ing the ray trajectory and polarizati o n stat e evolutions have 
been given in details in I Wang et al.l (|2014h . 

In Fig. [31 we show pulsar polarization emission beams 
for the X-mode, O-mode, total intensity, and the linear po¬ 
larization at six frequencies of v =0.8GHz, 1.0GHz, 1.4GHz, 
2.0GHz, 2.4GHz and 3.0GHz. The beams for low frequency 
emission are generally larger than those for high frequencies 
due to a larger opening angle from a higher altitude. More¬ 
over, the X-mode and O-mode dominate different parts of 
the pulsar beam, i.e., the X-mode dominates at two sides of 
the beam in £ direction while the O-mode dominates in (f> 
direction, as shown in the panels for Ix and Io in Fig. [3j 
Further, emission regions for the X mode and O mode are 
very distinct for lower frequency emission coming from a 
higher altitude. With increase of observing frequency, the 
distribution regions of the two modes gradually overlap, and 
hence serious depolarization happens, as shown by Ix — Io 
and fractional linear polarization of L% in Fig. [3] The dif¬ 
ferences are caused mainly by different corotation velocities 
of relativistic particles. When the corotation of relativistic 
particles is not considered, the distribut ion regions for tw o 
modes are almost completely overlapped (|Wang et al.ll2014h . 
The O-mode distribution regions are generally broader than 
those for the X-mode for the beams at higher frequencies, 
as shown in Fig. [3l which is caused by the dense plasma in 
deeper magnetosphere where serious refraction effect results 
in a prominent outward shift of the O mode beam. 

When a given sight line cuts across a pulsar beam, it 
will find polarized pulse profiles at a series of frequencies, as 
shown in Fig.[4l For the central cut at £ = 30°, the observed 
pulse profiles are predominantly of the O mode except for 
a small central phase region for the X mode. Emission at 
a lower frequency, like v — 0.8GHz, is highly linearly po¬ 
larized. With the increase of frequency, the fractional linear 
polarization, L%, gradually decreases and the profile width, 
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1 2 

y(GHz) 

Figure 5. Evolution of fractional linear polarization L% and pro¬ 
file width Wio with frequency. The fractional linear polarizations 
are calculated for the profiles in Fig. [4] with the sight lines of 
C = 30°, 31°, 32° and 33° at frequencies v = 0.8, 1.0, 1.4, 2.0, 2.4 
and 3.0 GHz. The red solid line represent the fractional linear po¬ 
larizations for the entire beams in Fig. [3] with the total intensity 
above 1% the peak intensity. The model parameters are the same 
as those in Fig. [3] i.e., a = 30°, 7 = 500, 77 = 100, B * = 10 12 Gs 
and P = Is. 


Wio, narrows down correspondingly, as shown in Fig. [5] As 
observing frequency further increases to about 3.0GHz, the 
pulse profile keeps narrowing down, but the fractional lin¬ 
ear polarization increases. This is because the refraction ef¬ 
fect becomes stronger as the plasma density increases, which 
dominates over the influence of rotation in determining pul¬ 
sar linear polarization. The frequency dependencies of frac¬ 
tional linear polarization and profile width are similar for 
different sight lines, such as £ = 31°, 32° and £ = 33°. More¬ 
over, the fractional linear polarization of the entire pulse 
beam exhibits a similar frequency dependence as those for 
the above discussed pulse profiles, which is outlined by the 
red solid line in Fig. [5] 

The frequency dependencies of fractional linear polar¬ 
ization and profile width are also calculated for plasmas 
with various Lorentz factors 7 and multiplicity parameters 
77 , as shown in Fig. [ 6 ] It can be seen from Fig. [ 6 ja) that 
profile widths decrease with the increase of frequency, but 
the fractional linear polarization decreases first and then 
may increase, similar to those in Fig. 0 except that the 
lower frequency emission is generated by relativistic parti¬ 
cles with a smaller Lorentz factor. If particles have larger 
Lorentz factors, emission of a given frequency will be gen¬ 
erated in higher magnetosphere. Because the rotation has 
already separated the distribution regions for the X and O 
modes above a certain height, emission above the height will 
be highly polarized with an almost constant polarization de¬ 
gree, as shown in Fig. |5{b). We can know from Fig. EUc) 
and (d) that high frequency emission generated in a dense 
plasma with a large 77 in deeper magnetosphere tends to 
have a larger fraction of linear polarization than those in 
less dense plasma with a small r\. This is caused by mode 
separation induce d by the O-mode refra ction, the effect is 
proportional to 77 (Bar nard h Aronslll986h . 

In summary, with curvature radiation mechanism, the 
radiation keeps highly linearly polarized at low frequen¬ 


ce 

j ^ 



0.5 1 

z/(GHz) 

(a) 7 = 400 , 77 = 100 


1 2 
v(GHz) 

(b) 7 = 600 , 77 = 100 




z/(GHz) z/(GHz) 

(c) 7 = 500, 77 = 10 (d) 7 = 500, 77 = 1000 

Figure 6. Same as Fig. [5] except for particles with different 
Lorentz factor 7 or density 77. 


cies. The polarization gradually decreases with observing 
frequency, but may increase at a very high frequency. The 
rotation plays an important role in keeping high degree of 
linear polarization, and can cause the decrease of fractional 
linear polarization in higher magnetosphere. Serious refrac¬ 
tion in a deep magnetosphere, which is significant when the 
plasma is dense and particles have a small Lorentz factor, 
will lead a larger fraction of linear polarization for high fre¬ 
quency emission. 


3.2 Emission generated from a region of the same 
height 

We have investigated the frequency dependence of pulsar lin¬ 
ear polarization emitted at the characteristic frequency z'cr- 
Curvature radiation from a relativistic particle or a particle 
bunch has a spectrum which is a power law with an index 
of 1/3 at the lowest frequency and decreases exponentially 
at high frequencies, not just at the frequency of ^cr- Hence, 
emission of various frequencies may be produced from the 
region of the same height in a pulsar magnetosphere. 

For investigations conducted here, it is assumed that 
emission of different frequencies is generated from the re¬ 
gion of the same height as the emission of ^cr = 1.0GHz. 
For a given i/c r, the shape of curvature radiation spectrum 
can be uniquely determined. Therefore, emission of differ¬ 
ent frequencies has the same beam pattern. However, due 
to the influence of frequency dependent propagation effects, 
the final beam patterns will be modified. As done in the 
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Figure 7. Same as Fig.[3]for polarization beams but with different parameters. Panels (a) are for the plasma with 7 = 500 and 77 = 100. 
Emission at frequencies of 0.4 and 6.0GHz is generated from the same height region as the emission of v cr = 1.0 GHz. Panels (b) are for 
the plasma with 7 = 400 and 77 = 1000. The other parameters used for model calculations are the same as those in Fig. [3] i.e., a = 30°, 
B* = 10 12 Gs and P = Is. 
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z^(GHz) 

(b) 7 = 400, 77 = 1000 


Figure 8 . Same as Fig. [5] but for (a) Lorentz factor 7 = 500 
and plasma density 77 = 100, and (b) Lorentz factor 7 = 400 and 
plasma density 77 = 1000. The fractional linear polarization and 
profile width are calculated by cutting the beams in Fig. [7] 


previous section, we calculate intensity distributions of the 
two mode within the entire pulsar beam, and then polarized 
pulse profiles for various frequencies and sight line angles, 
and finally find the corresponding frequency dependencies 
of fractional linear polarization and profile width, as shown 
in Figs. H and 0 

For moderate plasma conditions with a Lorentz factor 
of 7 = 500 and a density of 77 = 100, the wave mode dis¬ 
tributions within a pulsar beam are similar for various fre¬ 
quencies, as shown in Fig. [7| (a). The profile widths and 
fractional linear polarizations almost do not change with 
frequency, as shown in Fig. \8\ (a). This kind of frequency 
dependence is caused by the fact that the effects of rotation 
are the same for emission of different frequencies, because 
they are generated from the same height region, and the O 
mode components are less affected by the refraction effect. If 
the plasma has a smaller Lorentz factor of 7 = 400 but with 
a larger density of 77 = 1000 , the distribution regions for the 
X and O modes are overlapped, and the refraction effect be¬ 
comes much more significant, especially for lower frequency 
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Figure 9. Same as Fig. [3] except for particles with a gaussian 
distribution of energy, with 7 m = 500, and <r 7 = 50. 

emission ([Barnard &; Aronsl Il986l : I Wang et~aD l2014h . It re¬ 
sults in broader distribution regions for O mode than those 
for X mode, as shown in Fig. [7] (b). Compared with emis¬ 
sion at 6.0GHz, the O mode distribution region for emis¬ 
sion at 0.4 GHz is even much broadened. Therefore, both 
the profile width and fractional linear polarization decrease 
with frequency. This frequency dependence is due to that 
polarization features at higher frequencies are generated in 
deeper magnetosphere within a small altit ude range. O ur re¬ 
sult is contrary to the argument given bv iMcKinnonl (1 997 ) 
that the refraction causes average pulse profiles to depolarize 
while the profile width remains a constant. 


3.3 Emission generated by particles with 7 
distribution 

The frequency dependence of pulsar emission has been con¬ 
sidered for curvature radiation of relativistic particles with 
a single 7 in the regions with either a range of heights or 
the same height. Relativistic particles within a pulsar mag¬ 
netosphere generated by th e spa rking process may have an 
energy distribution ([Medin &; Lail I201Q ) . which we assume 
to be N(y) ~ exp [—(7 — y rn )/2a? f \. Doing the same as pre¬ 
vious subsections, we get intensities and polarizations of 
wave modes within the entire pulsar beam for each discrete 
Lorentz factor, and then sum them incoherently for total 
radiation. The polarized pulse profiles and their frequency 
dependencies are calculated, as shown in Fig. [9] which are 
similar to those for the relativistic particles with a single 7 
shown in Fig. [3j The frequency dependence here becomes 
less steep. Depolarization, which is caused by addition of 
emission by relativistic particles with different Lorentz fac¬ 
tors, leads emission at 3.0GHz to be less polarized compared 
to emission produced by particles with a single 7 . 


4 COMPARISON WITH OBSERVATIONS 

We have numerically simulated wave mode distribution re¬ 
gions and polarization beams, calculated polarized pulse 
profiles and their evolution with frequency for the conal- 
shaped density model with various plasma conditions. We 
have demonstrated the frequency dependent depolarization 
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Figure 10. The profile widths, VF 10 , and fractional linear polar¬ 
izations, L%, at a series of frequencies. The data are collected 
from literatures as listed in Table [l] The dotted lines represent 
model fitting. 


without lose of generality. In general, the fractional linear 
polarization together with profile width decrease with fre¬ 
quency. In order to check the validity of our simulations, it 
is worthwhile to compare our results with pulsar polariza¬ 
tion observations. 

Fig. [10] shows the changes of observed profile widths 
and fractional linear polarizations for six pulsars at a series 
of frequencies, with the data collected in table [l] in appendix 
from litera tures. These pulsars are of different mo rpholog- 
ical types ([RankinlIl983l : iLvne &; Manchester! Il988h . PSR 
J0528+2200 and J1136+1551 are conal double pulsars, 
each of which has two well separated conal components and 
smoothly varied position angle curves. PSR J0543+2329 
and J2022+5154 belong to the type of partial cone, which 
has only one conal component and represents one side of 
conal emission as indicated by position angle curve. PSR 
J2048-1616 is a triple-component pulsar. It has two conal 
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components at the outer parts of pulse profiles and one core 
component in middle. With increasing frequency, the cen¬ 
tral core componentgradually disappears as it has a steeper 
spectrum riRankinllTflSah . PSR J0452-1759 has multiple 
components. 

The evolution of profile width and fractional linear po¬ 
larization with frequency can be interpreted for these pul¬ 
sars by combing the curvature radiation processes together 
with propagation effects within pulsar magnetosphere. The 
widths of all these pulsars decrease with frequency. A con¬ 
stant plus a power law, Wio = W c + (yl v o) a , ca n accou nt 
for the frequency dependence, as pointed out by iThorsettl 
f 199 li b The frequency dependence of pulse width or beam 
width can be understood in the frame of curvature radi¬ 
ation mechanism by involving the density and energy dis¬ 
tributions o f relativistic particles within a pulsar magneto¬ 
sphere (I Wang et al .1120131 ) . The fractional linear polarization 
is found to decreases with frequency as well. A constant mi¬ 
nus a power law, L — L c — (v/iso) 1 *, can be used to fit the fre¬ 
quency dependence, as indicated by dotted lines in Fig. 1101 
which could result from depolarization by addition of the 
X mode and O mode as discussed in section 3.1 and 3.2. 
However, the fractional linear polarization of PSR J0452- 
1759 decreases first and then increases, and the increase of 
fractional linear polarization may be due to the significant 
refraction effect, as discussed in section 3.1. 


5 DISCUSSIONS AND CONCLUSIONS 

In this paper, we have investigated the frequency depen¬ 
dence of pulsar linear polarization by considering the cur¬ 
vature radiation process together with propagation effects. 
The X and O mode distributions within the entire pulsar 
beam, polarized pulse profiles and their evolution with fre¬ 
quency are calculated for emission at different frequencies 
originated from regions of different or the same height in a 
pulsar magnetosphere. We find the following conclusions: 

(i) Rotation can lead to the separation of distribution 
regions for the X and O mode within a pulsar beam, the 
effect of which is much more significant for low frequency 
emission generated at higher magnetosphere. 

(ii) The fractional linear polarization tends to approach 
a constant degree towards low frequency limit, as the distri¬ 
bution regions of the two modes have almost been separated 
above certain heights. 

(iii) The fractional linear polarization generally de¬ 
creases with frequency. Because the higher frequency emis¬ 
sion is generated from a lower altitude, where rotation in¬ 
duced mode separation is not significant. 

(iv) The fractional linear polarization could increase at 
higher frequency due to the significant O mode refraction. 

(v) If emission originates from almost the same height 
region, the O mode refraction will lead to the decrease of 
profile width and fractional linear polarization with fre¬ 
quency. The refraction influence is much more significant 
if the plasma has a large density or a small Lorentz factor. 

(vi) Depolarization caused by the randomization of po¬ 
larized emission from different heights of pulsar magneto¬ 
sphere is obvious for high frequency emission generated from 
the hot relativistic plasma. However, it does not change the 


basic tendency of which the fractional linear polarization 
decrease with frequency. 

The model simulations in this paper may help us to 
understand some statistical features of pulsar linear polar¬ 
ization. For short period pulsars, relativistic particles within 
the magnetosphere may suffer much more significant rota¬ 
tion influences. Hence, the X mode and O mode distribution 
regions tend to be separated and the pulsars may have higher 
linear polarization. The relativistic particles of high E pul¬ 
sars may be accelerated to a very higher energy (larger 7 ). 
Their emission can be generated from a much higher mag¬ 
netosphere, where the significant rotation will also lead to a 
large fraction of linear polarization. 

In our calculations, the frequency dependence of emis¬ 
sion region is determined by referring to the characteristic 
frequency of curvature radiation, and emission intensities are 
closely related to local plasma density. However, emission in¬ 
tensities depend on the detailed coherent process within pul¬ 
sar magnetosphere, which is not yet resolved now. Moreover, 
our calculations simply assume that streams of relativistic 
plasmas are cold and distributed in a conal shaped region, 
which travel along the curved dipole magnetic field lines 
and corotate with a pulsar magnetosphere. However, the ac¬ 
tual energy and density distributions of particles are un¬ 
known, which can be hot and with irregular density shapes 
([Medin fc Lail 12010 1. Furthermore, due to the influences of 
rotation and polar cap current, pulsar m agnetosp here may 
be of distorted dipole form (jKumar &; Gangadharal l2013h . 
which may affect pulsar polarization states and is not yet 
considered in this paper. 

Pulsars exhibit rich diversities of polarization profiles, 
some of which are not consistent with our model predictions. 
For example, the profile width of PSR J1239+2453 decreases 
with observing frequency, but its fractional linear polariza¬ 
tion increase first and then decrease; the profile widths of 
PSR J0601-0527 decrease but their fractional linear polar¬ 
izations increase; the profile width and fractional linear po¬ 
larization for PSR J2113+4644 increase. Moreover, individ¬ 
ual profile component s of a given pulsar ca n behave differ¬ 
ently with frequency (1 Johnston et al.l [2QQ8h . For example, 
the highly polarized profile component of PSR J0922+0638 
remain highly polarized at high frequencies. These diverse 
polarization features need further investigations. 
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Table 1. The profile width, W\o, and the fractional linear polar¬ 
ization, L%, at various frequencies for six pulsars. 


PSR Jname 

Freq. (GHz) 

WioO 

L% 

Ref 

J0452-1759t 

0.23 

21.1 

27 

1 


0.4 

19 

23 

1 


0.6 

18.5 

17 

1 


0.631 

18 

10 ± 2 

6 


0.92 

18 

12 

1 


0.95 

18 

13 ± 1 

10 


1.4 

18.3 

12 

1 


1.6 

18.5 

13 

1 


1.612 

18 

15 ± 1 

7 


1.72 

- 

16.9 ± 3 

8 


2.65 

- 

13.4 ± 3 

8 


4.85 

- 

19 ± 2 

11 

J0528+2200 

0.23 

22.6 

36 

1 


0.4 

20.8 

42 

1 


0.408 

- 

50 

4 


0.6 

20.4 

35 

1 


0.61 

- 

60 

4 


0.92 

19.4 

34 

1 


1.4 

18.8 

27 

1 


1.41 

- 

40.4 ± 0.2 

12 


1.6 

18.8 

27 

1 


1.71 

17.9 ± 0.7 

30.9 ± 1.5 

12 


1.72 

- 

33.5 ± 3 

8 


2.65 

- 

23.9 ± 3 

8 


4.85 

- 

20.6 ± 0.6 

12 

J0543+2329 

0.23 

32.0 

73 

1 


0.4 

27.6 

75 

1 


0.43 

- 

79 

9 


0.6 

24.3 

64 

1 


0.92 

24.2 

52 

1 


1.4 

25.7 

43 

1 


1.41 

- 

46.9 ± 0.9 

12 


1.5 

23.3 

43.2 ± 1.0 

13 


1.6 

26.5 

36 

1 


1.72 

- 

39.2 ± 3 

8 


2.65 

- 

27.4 ± 5 

8 


4.85 

- 

29.7 ± 1.2 

12 


8.75 

- 

17.3 ± 7 

8 


10.45 

22.5 ± 4.2 

6.2 ± 7.1 

12 


10.55 

- 

15 ±2 

14 

J1136+1551 

0.151 

_ 

35 

4 


0.23 

13.9 

25 

1 


0.24 

- 

35 

4 


0.4 

13.0 

23 

1 


0.408 

- 

30 

4 


0.41 

- 

33 

5 


0.6 

12.1 

20 

1 


0.61 

- 

25 

4 


0.638 

12 

24 ± 2 

6 


0.92 

11.3 

16 

1 


0.95 

11 

13 ± 4 

10 


1.4 

11.4 

9 

1 


1.41 

- 

15.8 ± 0.4 

12 


1.6 

10.5 

16 

1 


1.612 

10.5 

12 ± 1 

7 


1.665 

- 

12 

5 


1.71 

- 

12.1 ± 0.2 

12 


1.72 

- 

12.6 ± 4 

8 


2.65 

- 

8.4 ± 3 

8 


4.85 

- 

10.4 db 0.3 

12 


8.4 

8 

4 db 4 

3 


8.75 

- 

2.8 db 3 

8 


10.45 

- 

7.8 db 23.9 

12 


10.55 

- 

< 5 

14 


0.151 

- 

25 

4 

0.24 

- 

40 

4 

0.4 

18.6 

27 

1 

0.4 

19 

35 ± 2 

2 

0.408 

- 

35 

4 

0.41 

- 

40 

5 

0.6 

17.4 

21 

1 

0.61 

- 

55 

4 

0.638 

16 

37 db 2 

6 

0.92 

16.4 

24 

1 

1.4 

16.1 

21 

1 

1.5 

15.1 

37.5 db 0.1 

13 

1.6 

15.5 

18 

1 



1.612 

1.665 

1.72 

2.65 

4.85 

10.55 

16 

14.2 db 0.4 

31 ± 2 

30 

39.7 ± 3 
21.3 db 5 

9.2 db 1.6 

<C 5 

7 

5 

8 

8 

12 

14 

J2022+5154 

0.23 

24.8 

49 

1 


0.392 

- 

50 

5 


0.4 

20.5 

65 

1 


0.6 

19.7 

45 

1 


0.92 

17.8 

45 

1 


1.4 

16.9 

40 

1 


1.41 

- 

28.6 ± 0.1 

12 


1.6 

17.4 

32 

1 


1.665 

- 

34 

5 


1.71 

- 

38.8 ± 0.3 

12 


1.72 

- 

39.2 ± 6 

8 


2.65 

- 

33 ± 3 

8 


4.85 

- 

22.7 ± 0.3 

12 


8.75 

- 

25.7 ± 4 

8 


10.45 

12.7 ± 2.6 

29.6 ± 6.3 

12 


10.5 

- 

25 ± 1 

14 


14.8 

- 

15.9 ± 6 

8 


Notes. f The widths of PSR J0452-1759 are exceptionally measured at 50% 
of the peak intensities. References for previous polarization observations: 


1 Gould &; Lvne (| 19981^), ^Hamilton et al. (19771), ^Johnston et al.l ( 

2006 

), 

' Lvne et al.l dl97l|), ^Manchester! (197ill'), ‘iMcCulloch et al.l ( 

1978 

), 

Manchester et al.l il980l), ^Morris et a 

4.1 (198l|), ^Rankin & Benson! 

198 lh , iL lvan Ommen et al. Jl997l), 11 

von Hoensbroech et al.l (1 

119981 


von Hoensbroech Xilouris! (|l997|), i;:i 

Weltevrede & Johnston! 1 

120081 

), 


J2048-1616 
































